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^As  part  of  an  effort  to  determine  electron  density  profiles  in  the  auroral 
region  of  the  ionosphere,  a  method  has  been  developed  for  directly  calculating 
real  height  from  virtual  height  data  obtained  at  high  latitudes.  The  essential 
ingredient  of  the  method  is  the  representation  of  the  curved  portion  of  the 
ordinary  ray  dispersion  curve  by  an  analytic  form  which  reduces  the  virtual 
height  -  real  height  integral  relation  to  a  form  of  Abel's  integral  equation. 
Numerical  checks  indicate  that  the  present  method  gives  an  accurate  represen¬ 
tation  of  the  dispersion  curve.  The  correctness  of  the  present  procedure  has  - 
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^been  verified  for  the  case  of  a  linear  electron  density  profile. 

The  present  method  will  be  applied  to  the  numerical  calculation  of  real 
height  electron  density  profiles  in  the  auroral  region  from  virtual  height  data. 
The  accuracy  and  convenience  of  the  present  method  will  be  assessed  by 
comparison  of  the  results  with  those  obtained  from  standard  programs. 
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A  Method  for  Direct  Determination  of  Real 
Height  From  Virtual  Height  Data  for  the 
Auroral  Region  of  the  Ionosphere 


1.  INTRODUCTION 

For  the  case  of  no  magnetic  field,  the  group  index  of  refraction  is  of  fairly 
simple  form  and  the  virtual  height-real  height  Integral  relation  may  be  expressed 
as  a  special  form  of  Abel's  integral  equation  and  solved  for  real  height  as  a  direct 
function  of  virtual  height.  1 

In  the  presence  of  a  magnetic  field,  the  group  index  of  refraction  becomes 

quite  complicated  and,  in  general,  the  group  height- real  height  integral  relation 

cannot  be  reduced  to  a  form  of  Abel's  integral  equation.  The  approach  generally 

2  3 

used  here  is  the  laminar  method  *  in  which  the  electron  density  profile  is  divided 
into  a  series  of  slabs  in  each  of  which  average  values  are  used  to  describe  the 
properties  of  the  slab.  The  integral  relation  then  reduces  to  a  set  of  simultaneous 
equations  whose  coefficients  may  be  calculated,  giving  the  real  height  as  a  function 

of  virtual  height.  Improvements  to  the  laminar  method  have  been  provided  by 

4  5 

Titheridge  by  use  of  a  single  polynomial  or  a  series  of  overlapping  polynomials. 

For  the  particular  case  of  longitudinal  propagation,  the  group  index  of  refrac¬ 
tion  becomes  fairly  simple  and  the  group  height- real  height  relation  may  be  reduced 
to  a  form  of  Abel's  integral  equation.  Since,  in  the  auroral  region,  the  angle  0 
between  the  magnetic  field  and  the  incident  ray  is  less  than  about  30*  for  vertical 


(Received  for  publication  13  November  1979) 

(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  19.) 


propagation,  it  is  natural  to  inquire  whether  the  group  height-real  height  relation 
may  be  expressed,  at  least  approximately,  in  terms  of  Abel's  integral  equation. 

For  the  extraordinary  ray.  Figure  1,  the  index  of  refraction  remains  virtually 
unchanged  for  $  less  than  about  30° 6  and  the  real  height  in  the  auroral  region  may 
be  determined  from  an  Abel  integral  equation  in  a  fairly  straightforward  manner. 


Figure  1.  Sketch  of  Ordinary  and  Extraordinary 
Dispersion  Curves  for  0  Less  than  30* 


For  the  ordinary  ray,  however,  as  seen  in  Figure  1,  the  situation  is  somewhat 
more  complicated.  Although  M2  is  almost  linear  over  a  portion  of  X  =  fN  / f2 , 
where  jjis  phase  index  of  refraction,  f^  is  plasma  frequency  and  f  is  probing  fre¬ 
quency,  the  curve  bends  rapidly  at  the  larger  values  of  X  to  give  a  reflection  level 
at  X  *  1.  It  is  found  that  the  nonlinear  portion  of  the  curve  may  be  fitted  by  a 
simple  power  law  for  which  an  Abel  integral  equation  obtains.  Thus,  the  dispersion 
relation  for  the  ordinary  ray  may  be  regarded  as  made  up  of  two  segments  for  each 
of  which  the  real  height  may  be  obtained  directly. 


6.  Davies,  Kenneth  (1989)  Ionospheric  Radio  Waves.  Blaisdell  Publishing  Company, 
Waltham,  Massachusetts,  Chapter  b,  pp  136-135. 


2.  ANALYSIS 


2.1  Extraordinary  Ray 

For  the  extraordinary  ray  the  index  of  refraction  in  longitudinal  propagation 
la  given  by 

M2-l-py  (1) 

where 

fH 

Y  *  -j-  is  the  gyromagnetic  ratio 

and  f„  is  the  gyromagnetic  frequency.  As  indicated  previously,  the  dispersion 

Cl 

relation  given  by  Eq.  (1)  remains  virtually  unchanged  for  0  less  than  about  30*. 

We  may  therefore  use  Eq.  (1)  in  the  analysis  of  the  relation  between  group  height 
and  real  height  for  the  auroral  region. 

The  virtual  height  h  '  (f ),  which  is  a  function  of  input  frequency  f,  is  related  to 
the  real  height  h  by  the  integral 

hr 

h'(f)=  J  p'dh  (2) 


where  h  is  real  height,  hr  is  the  height  at  reflection,  that  is,  when  4=  0,  and  4' 
is  the  group  index  of  refraction  and  given  by 

4  =  gy-<4f)=4+f  .  (3) 


The  group  index  may  now  be  expressed  explicitly  as  a  function  of  X  and  Y.  For  our 
purposes,  however,  it  will  be  more  convenient  to  use  the  first  form  of  Eq.  (3).  We 
may  then  write  Eq.  (2)  as 


/  ^  (4f)  dh  =  h,(f)  * 


(4) 


Noting  that  at  the  upper  limit  hr  the  index  4  vanishes,  we  can  take  the  derivative 
outside  the  integral  sign  and  write 


d 

3T 


r 

l 


lit  dh 


h'(f)  . 


(5) 
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We  may  now  integrate  with  respect  to  f  to  obtain 


r  f 

/  f(l  -  r\)il2  dh  =  /  h'(f)df 
0  "  f. 


H 


(6) 


where  the  lower  limit  of  f  is  set  at  fjj  since  X  is  negative  in  the  range  0  to  f  when 
Y  >  1.  Noting  that  Y  is  independent  of  the  variable  dh,  we  may  transpose  the  term 
in  1  -  Y  to  the  right-hand  side  and  write  Eq.  (6)  in  the  form 


(f2  -ffH  -f2)172  dh*  H'(f)  . 

(7) 

f 

H'(f)  =  (1 -Y)172  J  h'(f)df. 

(8) 

To  reduce  Eq.  (7)  to  the  form  occurring  in  Abel's  integral  equation  we  introduce 
new  variable  F  and  FN  where 


F  =  f2  -  ffH  . 

F  =  f  2 
FN  *N  * 


which  allows  us  to  cast  Eq.  (7)  in  the  form 


^(F  -  $  )1/2 


u  d£  =  H'(f) 


(9) 


where  {  =  FM  is  a  running  variable  with  maximum  value  F  corresponding  to  the 
reflection  frequency  and  u  *  . 

We  may  express  Eq.  (9)  in  the  canonical  form  of  Abel’s  integral  equation  by 
differentiating  both  sides  of  Eq.  (9)  with  respect  to  F.  However,  the  solution  can 
also  be  obtained  by  direct  use  of  the  Laplace  transform,  indicating  at  the  same  time 
the  simplicity  and  utility  of  the  Laplace  transform  technique  in  this  type  of  integral 
equation.  Applying  the  Laplace  transform  to  Eq.  (9)  we  obtain 


V7  1 


T~  ^TT2 


u(s)  *  H<s) 


(10) 


where  the  bar  indicates  the  transformed  variable,  and  s  the  transform  parameter, 
that  is. 


8 


and  we  have  made  use  of  the  relation  for  the  Laplace  transform  L  of  a  convolution 
integral  for  two  functions  f({)  and  g({ )  in  terms  of  the  individual  transforms 


fls)  g(s)  *  L  J  f(£  -  ij)  g(q)  d>)  . 
0 


(12) 


Solving  Eq.  (10)  for  u(s)  we  obtain 
u(s)  =  —  s3/2  Ifts)  . 

V* 


(13) 


Since  the  inverse  Laplace  transform  for  s  does  not  exist  for  n  >'0  except  when  n 
is  an  integer  we  write  Eq.  (13)  in  the  form 


u(s)  =  — —  —  i  ra1  ff(i) 

VF  a17- 


(14) 


whose  inverse  transform  is  readily  calculated  as 


*  dl*  0  (5  -F)1/Z 


(15) 


to  the  plasma  frequency  FM  at  reflection.  Since  u  =  ^  we  can  write  Eq.  (15)  as 


where  F  is  now  the  running  variable  and  £  is  the  maximum  value  of  F  corresponding 

dh 

cav  x  uiULc  u  —  ^ 

.  2d  f  H  (F) 

h(*  >  ~  ?  ar  1  ;r  dF 


(16) 


0  (£  -  F) 


which  can  be  further  written  in  the  forms 
F 

h(F  )  =  |  /  - 1 — rjy  d?p(F>  d  F  . 

v  ’  0  (Fy-F)172  ** 


x 

x  2  rV  1  dH'(f) 

Mfv)-?-l  tl  — df- 
*H 


(17) 


(18) 


where  Fv  or  fy  is  the  maximum  value  of  F  or  f  corresponding  to  the  reflection 
height,  and 


r1  ■  K2  ■  'v'h  -  «2  -  «H>  • 


Carrying  out  the  differentiation  we  obtain 


*v=f  / 


(1-Y) 


1/2 


h'  (f)  df  +  i  / 


A  T^Y  H'(f)df  • 


(19) 


‘H 


H 


2.2  Ordinary  Ray— Linear  Segment 

As  indicated  previously,  the  ordinary  ray  remains  almost  linear  over  some 
range  of  X  and  then  begins  to  bend  rapidly  toward  p  =  0  at  X  =  1.  We  shall  there¬ 
fore  consider  the  ordinary  ray  as  made  up  of  a  linear  and  a  curved  segment  joined 

o 

at  the  common  point  X  =  X j.  The  point  X ^  is  chosen  at  the  point  at  which  the  fi 
curve  starts  to  bend  away  noticeably  from  the  straight  line  portion  of  the  curve. 

The  contribution  of  the  straight  line  portion  of  the  curve  may  be  obtained  in 
the  same  manner  as  for  the  extraordinary  ray.  However,  to  use  the  Abel  integral 
relation  it  is  necessary  that  the  upper  limit  refer  to  the  point  X  at  which  n  vanishes. 

n 

The  lower  limit  may  be  arbitrary,  so  long  as  H'  (f)  vanishes  at  the  lower  limit, 
a  condition  which  is  automatically  satisfied.  We  therefore  extend  the  straight  line 
to  H  =  0  at  X  =  1+Y  and  calculate  the  contribution  to  the  real  height  of  the  straight 
line  portion  from  0  to  Xj  as  the  difference  of  the  contributions  from  0  to  1  +  Y  and 
from  Xj  to  1  +  Y. 

For  the  straight  line  portion  between  0  and  1  +  Y  the  real  height-virtual  height 
relation  is 


* 

f  (|if)dh=  /  h'(f)df 

0  o 


(20) 


.  .  2  . 

but  H  is  now  given  by 


4=1- 


T+T 


(21) 


The  lower  frequency  limit  can  here  be  taken  as  zero  since  X  remains  positive  for 
the  ordinary  ray  when  Y  becomes  greater  than  unity.  Following  the  same  procedure 
as  for  the  extraordinary  ray  we  obtain 
h 


/ 


(f2+ffH*fN2) 


1/2 


dh  =  H'(f)  . 


(22) 


which  can  be  written  as 


/  (F  -  4  )1/2  udi  =  H'(F)  . 
0 


(24) 


where 


f  =  r  +  «H , 


«  =FN=fN  ' 


The  real  height  is  then  given  by 
Fv 

u  ,  _  2  r  1  dH'(F) 

0(Fv)  '  r  JQ  (F  _F)i/2  -3T  dF 


(25) 


or 


VV  ■  I  /  i  T8 


(26) 


where 


X2  ■  'v2  *  fv  <H  -  «2  +  "h1 


Carrying  out  the  differentiation  we  obtain 


2  r v  (l  +  Y)1^2 


w  =  !/  h'(ndf-l/  A  TtT  «'(»«■ 


(27) 


For  the  interval  between  Xj  and  1  + Y  the  height  relation  is 
h 


r  2  2  1/2 


/  (r+ffjj-fj/)  dh  =  H'  (f)  , 
h! 


(28) 


f 


h:  <f)  =  d+Y)1/2  /  h'(f)  df 
1 


(29) 


11 


fcgBB SB 


J 


where  the  subscript  1  corresponds  to  Xj  and  the  frequency  f  j  is  given  by 


t\  +  flfH"  fNt 


(30) 


We  now  write  Eq.  (28)  in  the  form 
F 

/  (F  -  i  )1/2  ud*  =  H  '  (F) 

which  may  be  solved  for  h^  (M  to  yield 
F 


h1(Fv) =  I 


l 


(fv-f) 


.  d  H. '  (F) 

<n?  dF 


(31) 


(32) 


or 


(f .  2 ;  i  dHi'<f>  df 

h  (fv)  "  J  1  R  cTT -  df  ‘ 


(33) 


The  height  hQ1  between  0  and  X1  is  accordingly  given  by 


h01  =  h0  ‘  hl 


(34) 


2.3  Ordinary  Ray— Determination  of  Curved  Segment 

We  represent  the  curved  portion  of  the  dispersion  curve  in  the  interval  Xj  to  1 


by 


U2  =  a<l-X)“  H2 


(35) 


where  a  and  a  are  parameters  to  be  determined  and  4  j  is  the  index  of  refraction 
given  by  the  straight  line  portion  of  the  curve  at  the  point  Xj,  that  is, 

X, 


11  \  =  1  'TTY 


(36) 


The  value  of  a  is  between  0  and  1  and  approaches  zero  as  6  becomes  small. 
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2 

Since  (1  -X)  goes  to  zero  at  X=  1,  the  derivative  of  with  respect  to  Y 

vanishes  at  X=  1.  The  exact  dispersion  curve  is  also  known  to  be  independent 

2 

of  Y  at  X=  1,  so  that  the  form  for  fi  in  Eq.  (35)  has  the  proper  qualitative  depen- 

2 

dence  upon  Y.  In  fact,  a  computational  check  of  n  from  Eq.  (35)  as  a  function 

of  Y  was  in  good  agreement  with  the  exact  dispersion  equation, 
o 

The  slope  of  the  /i  -X  curve  given  by  Eq.  (35)  is  infinite  at  X=  1#  whereas  the 

g 

exact  curve  has  the  form,  near  X  =  1, 

H2  = - 2 —  <1-X)  (37) 

sin  9 

with  a  slope  of  magnitude - ij —  .  However,  since  9  is  less  than  30°,  the  slope 

sin  9  2 

from  Eq.  (37)  is  generally  large.  Thus,  the  u*  curve  given  by  Eq.  (35)  should 
generally  be  in  good  agreement  with  the  exact  dispersion  curve. 

2 

The  values  of  a  and  a  were  determined  by  fitting  the  value  of  u  in  Eq.  (35) 
to  those  of  the  exact  dispersion  curve  at  Xj  and  a  point  X^  between  Xj  and  1.  The 
results  obtained  for  9  =  10°,  20°  and  30°,  and  for  Y  =  0.  3  are  shown  in  Figure  2 
where  it  is  seen  that  the  fitted  curves  are  in  good  agreement  with  the  exact  curves. 
The  accuracy  of  the  fit  was  found  to  be  fairly  insensitive  to  a  moderate  change  in 
the  point  X2«  say  from  0.  9  to  0.  8.  In  addition,  although  the  slope  at  Xj  was  not 
determined,  the  value  from  Eq.  (35)  was  in  fairly  good  agreement  with  that  ob¬ 
tained  from  the  exact  dispersion  curve.  It  thus  appears  that  the  form  chosen  for 
2 

H  in  Eq.  (35)  is  a  useful  and  accurate  representation  of  the  curved  portion  of  the 
exact  dispersion  curve. 


2.4  Calculation  of  Real  Height  for  Curved  Segment  of  Dispersion  Curve 

2 

The  apparent  height  h'(f)  for  the  interval  X^  to  1,  with  fi  now  given  by  Eq.  (35) 
is  then 


i, 


fa^Mx(l 


a 

X)1  dh  = 


l 


h'  (f)  df 


(38) 


which  can  be  expressed  as 

•2  /-  9  9 >7  M 


r  a 

{ 1,2 -‘*v 


1 


a~- T 


dh  = 


=  (h' 


(f)  df  . 


(39) 


8.  Davies,  Kenneth  (1969)  Ionospheric  Radio  Waves,  Blaisdell  Publishing  Company, 
Waltham,  Massachusetts,  Chapter  5,  p  125. 
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The  terms  in  and  f  are  independent  of  real  height  h  and  we  can  accordingly 
write  Eq.  (39)  in  the  form 


(41) 


H 


jot- 1  . 

'(f)  =  -srn -  /  h'(f)df 

»W  fl 


Changing  to  the  variables  F  and  {  we  write  Eq.  (40)  in  the  canonical  form 
F  a 

’  (F) 


|  <F-{  )  ud*  =  H'l 


(42) 


which  yields  under  Laplace  transformation 


u(s>  r  (?  +  1) 


=  H(s)  . 


+  1 


(43) 


Solving  Eq.  (43)  for  u(s)  we  have 

a 


u(s) 


1 


+  1 


r<|  ♦  i) 


H(s) 


(44) 


which  may  be  written  as 


u(s) 


1 


r  (f  +  i  ai  “  ? 


R(s). 


The  inverse  transform  of  Eq.  (45)  is  then 


dh 

<nr 


2  « 
a  L 


T(l+f  )  d4  2  Fj  r(l-f  ) 


H'(F)dF 


(45) 


(46) 


which  may  be  written  as 
F 

.  Jf  a  v 

sin-^  f  1  dH'(F) 

h>  ■  -if-  4,  st- 


or 


h,  = 


.  n  a  v 

8inT  r  1  dH'(f) 

Ta  ]  ~  2_f2)a/2  37 - 


(47) 


(48) 


15 


where  we  have  used 


wa 

ro+f  >  ru-?)  -  ~ 

4  4  sin 


ra 


(49) 


Carrying  out  the  differentiations  we  may  write  Eq.  (48)  a? 
*  a  t 


sin 


‘1 


~T 

f7 


1  X1 

M  T 


/. 


df 


fl  <f  2-f2) 


a7T 


{?-■  *-  h'  (f)  -  (-if^  H' 
a  <*1 


(f) 


-  7  -iw-  H/  (f) 

(1  +  vr  Mj 


(50) 


3.  VERIFICATION  OF  REAL  HEIGHT  FORMULA  FOR  A 
UNEAR  ELECTRON  DENSITY  PROFILE 


Since  the  real  height  formulas  derived  here  are  of  the  same  general  type,  we 
shall,  for  convenience,  work  with  the  equations  for  the  extraordinary  ray.  For  a 
linear  gradient 


fN  =ah 


<51> 


where  fN  and  h  are  referred  to  suitable  initial  values  which  may  be  taken  as  zero. 

The  apparent  height  h'(f)  is  then  given  by 
h 

r 

h'(f)  =  3T  {  titdh 


/  ,  2  \  1/2 

*'<«=£/  I  Pity  -jr  J  dV 


(52) 


where  the  upper  limit  corresponds  to  the  value  of  at  which  H  vanishes.  The 
integral  yield 


[»  ll-VK3]. 

The  quantities  H'(f)  and  d  jy  —  are  then  given  by 


(53) 
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H'(f)  -  (l-Y)7  /  h'(f)df 

(54) 

(55) 

which  can  be  expressed  in  terms  of  F  by 


H'(f)  »  £  (l-Y)3/2t3 


dH'(f) 

ar? 


(i-vr 


f(2f  -  fH) 


dH'(F)  F2 

~ar  *  "a 


The  real  height  hr  is  then  obtained  from 


F 

2  fV  1  dH'(F)  Jf, 

r  ?  L  jtt  ~arr~  dF  * 


0  (Fy-F) 


F  1 

v  7 

k  2  r  F  2 

hr  fa  • 


0  (Fv-F) 


m 


dF 


Changing  the  variable  to  G*  F1^  yields 

Gv 

'  *  s„  757^7^  d0  ’ 


which  can  easily  be  evaluated  as 

h  -Is. 

hr  ~ST  ' 


(56) 


(57) 


which,  from  Eq.  (51),  is  identical  with  h^.  The  solutions  for  the  ordinary 
ray,  both  linear  and  curved  portions,  can  easily  be  verified  in  a  similar  manner. 


4.  SUMMARY  AND  CONCLUSIONS 

A  method  has  been  developed  for  directly  calculating  real  height  from  virtual 
height  data  in  the  auroral  region  of  the  ionosphere.  The  principle  approximation 


17 


I 


L 
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in  the  method  is  the  representation  of  the  curved  portion  of  the  dispersion  curve  for 
the  ordinary  ray  by  an  analytic  form  which  reduces  the  virtual  height-real  height 
integral  relation  to  a  form  of  Abel's  integral  equation.  Numerical  checks  show 
that  the  present  method  gives  an  accurate  representation  of  the  exact  dispersion 
curve.  It  is  therefore  felt  that  good  results  should  be  obtained  when  the  present 
procedure  is  applied  to  virtual  height  data.  The  correctness  of  the  real  height 
formula  has  been  verified  for  the  case  of  a  linear  electron  density  profile. 

The  present  method  will  be  applied  to  the  numerical  calculation  of  real  height 
from  virtual  height  data  and  its  accuracy  and  convenience  will  be  assessed  by 
comparison  of  the  results  with  those  obtained  from  standard  programs. 
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